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ABSTRACT 



Context. Photometry of Procyon obtained by the MOST satellite in 2004 has been searched for p modes by several groups, with sometimes 
contradictory interpretations. 

Aims. We explore two possible factors that complicate the analysis and may lead to erroneous reports of p modes in these data. 

Methods. Two methods are used to illustrate the role of subtle instrumental effects in the photometry: time-frequency analysis, and a search 

for regularly spaced peaks in a Fourier spectrum based on the echelle diagramme approach. 

Results. We find no convincing evidence of a p-mode signal in the MOST Procyon data. We can account for an apparent excess of power close 
to the p-mode frequency range and signs of structure in an echelle diagramme in terms of instrumental effects. 
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1. Introduction 

The search for acoustic oscillations in the star Procyon A 
(which we call hereafter simply Procyon) has been a key 
part of the young field of stellar seismolo gy for two decades , 
since the first claim of p-mode detection by Gel ly et al. (1986). 
In recent years, interes t grew rapidly with the announcement 
by iMartic et al] dl999h of global acoustic modes in Procyon 
detected in high-resoluti o n, tim e-resolved spectroscopy, con- 
firmed by iBarban et al.1 £l999). Interest turned to contro- 
versy with the non-detection of intensity oscillations by the 
MOST (Mic rovariability & Oscill ations of STars) space-borne 



photometer (Matthews et al.l |2004|) . whi le contempo r aneou s 



groundbased spectroscop ic c ampaigns bylMartic et al.l (12004|) . 
lEggenberger et al.l ( 2004 ) and lBouchy et al] d2004 ) appeared to 
confirm the presence of stellar oscillations. This non-detection 
triggered considerable sp eculation about non- stellar noise 
sources in the MOST data ( Christensen-Dalsgaard & Kjeldsenl 
20041: iBeddingeTa l. 2005) and independent analyses of those 



data claiming possib le evidence for p modes, but using different 
analysis methods (Regulo & Roca Cortesl 120051: Garcia et al. 
The contrast between the reported p-mode signal in lu- 
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minosity and integrated radial velocity also inspired theorists 
to refine models of c onvection and turbulence in Procyon (e.g., 
iRobinson etaDl2005l) . 

The MOST team has described the nature of the Fabry 
Imaging data on w hich the Procyon photometry is based 
dReegen et al.l l2006), including the recognised effects of stray 
light due to Earthshine modulated by the MOST satellite orbit 
and a beating effect of the science and startracker CCD timing 
in the early stages of the mission. 

We present an independent analysis of the MOST Procyon 
data compared to simulations, calling attention to how the tim- 
ing effects and orbital artifacts can be misinterpreted as evi- 
dence for p-mode oscillations. 

2. Time-frequency analysis 

The 32-day-long MOST photomet ric time series has a du ty cy- 
cle of 99% (for more details, see iMatthews et al.l 12004). This 
continuity of the data, currently unique to MOST observa- 
tions of stars and a key element of the motivation for future 
photometric space missions such as CoRoT and Kepler, has 
many important advantages. It eliminates the cycle/day aliases 
common to groundbased campaigns, even with multiple sites. 
Another benefit, which we exploit in this paper, is that these 
data make it possible to investigate variation of the frequency 
content of the signal as a function of time more sensitively than 
for other astronomical data sets. 

We have applied a time-frequency analysis, based on a Morlet 
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Fig. 1. Left: Time/frequency diagramme of Procyon data, showing excess power between 500 and 1200 //Hz during the first eight 
days of observation. Right: same analysis of the CCD temperature, showing the same excess of power 



wavelet (see iBaudin et al.1 Il994l) . to the time series of the 
MOST photometry and to time series of instrumental param- 
eters such as satellite pointing error and CCD temperature, 
which are included in the MOST science telemetry. The left- 
hand panel of Figure [T] shows the distribution of power in the 
MOST Procyon photometry as a function of both frequency 
and time. There is some excess power during the first 8 days 
of the observations, in the frequency range of about 500- 1200 
//Hz. The contrast in excess power is small and can be seen 
more clearly when the power is summed over frequency as a 
function of time, as shown in Figure |2] The average increase 
in the excess power is about 35%. The same time-frequency 
analysis was applied to time series of spacecraft and instru- 
ment parameters in the MOST science telemetry. Some (but 
not all) temperatures measured on board seem to show an ex- 
cess of power in the same ranges of time and frequency as in 
the stellar light curve. The time-frequency plot of the MOST 
science CCD temperature is shown in the right-hand panel of 
Figure [T] Other parameters, such as telescope-pointing wan- 
der, do not show such behaviour. However, this does not mean 
that on board temperature variations are the cause of the ob- 
served photometric power excess. MOST CCD temperatures 
(operating at around -35°C) are controlled to within 0.1°C and 
are monitored to an accuracy of 0.0 l°C The very small am- 
plitude of the temperature variations is insufficient to change 
the dark current or other characteristics of the CCD output to 
account for the excess seen in Figure [2 Furthermore, there is 
no one-to-one correlation in time and frequency. More likely, 
both the temperature and photometric variations are symptoms 
of a common cause. There is nothing in the MOST operations 
records or in other telemetry that can explain such a change on 
the seventh day of the observ ations. Amon g the d ifferent instru- 
mental effects described by iReegen et al.1 d2006|) . none seems 
to be clearly the cause of this observation. However, one must 
note that these seven days represent the first seven days of sci- 
ence operations of the MOST satellite. 

The smoothed Fourier spectrum of the first seven days of the 
Procyon photometry shows greater spectral power density up to 




Fig. 2. Time variation of the integrated power (from 400 to 
1600 //Hz), showing two different levels during the observa- 
tions. 

about 1200 //Hz than during the remaining 24 days (Figure [3]). 
This includes the frequency range where acoustic vibrations of 
the star are anticipated, both from theory and the spectroscopic 
observations. Therefore, this power contrast must be taken into 
account in any analysis looking for p modes, which should be 
sought preferentially in the second part of the photometric time 
series. 

3. Collapsed echelle diagramme analysis 

The acoustic spectrum of a star such as Procyon is expected to 
have modes with a roughly regular spacing in frequency (called 
the "large spacing" A ). One way to see this spacing (as well 
as fine spa cings betw een modes) is to construct an echelle di- 
agramme (IGredll98l[) . widely used in helio- and asteroseis- 
mology. This is done by cutting the Fourier spectrum into N c 
frequency intervals of width A, and then plotting each interval 
above the previous one, such that the x-axis is the frequency 
modulo A and the j-axis is the frequency. The regular spac- 
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Fig. 3. Smoothed Fourier spectra of Procyon data for the first 
seven days (black) and last 25 days (red) of observation 



ing of modes in the aymptotic regime yields clear, nearly ver- 
tical, ridges in the echelle diagramme. To enhance the signal- 
to-noise ratio (SNR) in an echelle diagramme, it is possible to 
collapse th e second dimens ion (the order of the echelle) as pro- 
posed by iKorzennikl (|2007|) . If the spacing A chosen is the ac- 
tual large spacing Ao of the p-mode eigenspectrum of the star, 
then a peak should appear or become clearer in this "collapsed 
echelle diagramme". To search for a p-mode spacing in data of 
low SNR, this can be repeated for a set of A values and a 2D 
collapsed echelle diagramme can be constructed. The result- 
ing diagramme, built with prewhitened spectra (by dividing the 
original spectra by the smoothed -200 jjHz window- spectra), 
is normalised as follows: 



Piy) = ( 



Piy) 

2N c S f 



(1) 



where s is a detection level defined below, P(v) is the nor- 
malized, collapsed echelle diagramme; p(v) is the original col- 
lapsed echelle diagramme having a x 2 with 2 N C S / degrees of 
freedom, N c being the number of frequency intervals collapsed; 
and S / is the smoothing factor used, if any. The normalisation 
factor (2N C S f) in Equation [Hallows one to compensate for the 
varying number of degrees of freedom in the collapsed echelle 
diagramme. 

We first set the a priori detection probability such that we 
can derive the detection level s using the cumulative proba- 
bility for the statistic s mentionned above. The detection level s 
( App ourchauxl2004l Equation 1) is then given by: 



P(s' >s,q) - 
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,q-l 



T(q) S« 



-u/S 



du 



(2) 



with q = N C S / and S is the mean noise level in the power spec- 
trum. The detection level ( s ) is set a priori to be 1%, i.e. there 
is a one percent chance that in a window of width N w bins, a 
peak due to noise will be higher than the level s. When there 
is no smoothing factor, the number of independent bins is N w . 
When a smoothing factor is applied, the number of independent 
bins decreases linearly with the smoothing factor. The level s 



for a single bin is then obtained by solving the following equa- 
tion: 



P(s' >s,q) = 0.0l-f 



(3) 



This approximation is valid because the probability of 1% is 
small compared to one. 

It should be noted that if a spacing Ao exists in the spectrum, a 
different spacing A Q defined as 



mAg = hAq, 



(4) 



with m and n being integers, will also be observed in an echelle 
diagramme. This artifact plagues any technique based on the 
detection of regularly- spaced peaks. 

We apply this method to simulated data in the next subsection, 
in the framework of pre paration for the CoRoT mission (see 
lAppourchaux et al.l 120061) in order to detect low-SNR modes 
with a regular frequency spacing. We then apply it to the MOST 
Procyon photometry. 

3.1. Simulated data 



The simulated data used here (see Baudi n et al. 2006) include 
an estimate of the granulation noise (the dominant component 
of the intrinsic stellar noise in the p-mode frequency range) 
and theoretically predicted p-mode amplitudes. The granula- 
tion noise is modelled in a simple way (a Lorentzian shape in 
the Fourier spectrum ) and the modes' amplitudes are based on 
Samadi et al.l (l2QQ3b . These amplitudes can be attenuated or en- 
hanced in order to produce a desired SNR. In the simulations 
here, we set the mode SNR (defined as (.Pmode + ^noise)/^noise) 
to be 1.3 . The input frequencies of the simulated modes cor- 
respond to a mean large spacing of Ao = 82yuHz, with a stan- 
dard deviation of 23jiHz. A smoothing factor correspond- 
ing to 3jiH z has been used. Figure H] shows the collapsed 
echelle diagramme for this simulation in which clear peaks 
above the normalised value of 1 .0 (indicating a 99% confidence 
level) appear for regularly- spaced frequencies having A ^ 41, 
82, 123, and 165 jjHz. These correspond to the input value of 
the large spacing A and the corresponding values A Q (defined 
in Equation|4]), demonstrating the strong potential of this ap- 
proach. 

3.2. MOST Procyon data 

The same type of collapsed echelle diagrammes (with the same 
smoothing factor) were computed for the MOST Procyon pho- 
tometry. They have also been applied to the window function 
of the MOST time series to assess its influence on the analysis. 
The gaps in the data represent only 7.5 hours in a total obser- 
vation of 32 days. This is illustrated by the window function 
(see Fig 13) in which the presence of very short gaps which oc- 
cur with some periodicity (related to the satellite orbit), induces 
very weak peaks. Despite their weakness, the regular spacing 
of these peaks in the Fourier spectrum can interfere with the 
search for regularly- spaced p modes in this spectrum. The left- 
hand panel of Figure [6] shows the collapsed echelle diagramme 
of the window function for the region of the Fourier spectrum 
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Fig. 4. Collapsed echelle diagramme for simulated data, show- 
ing the signature of regular spacing of ^ 82 jjHz and multiples 
of this spacing as defined in Equation |4] 
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Fig. 5. Power spectrum of the window function of the obser- 
vation, almost perfect but presenting some weak but regularly- 
spaced peaks 



where stellar p modes are expected. One can see clear verti- 
cal and oblique structures: the vertical ones are due to large 
peaks (whatever their origin) in the first frequency interval used 
to build the diagramme, and the oblique ones to similar large 
peaks, present in subsequent frequency intervals. Clear peaks 
above the 99% confidence level appear at ^165, 110, 83, and 
55 //Hz. The first of these frequencies is the orbital frequency 
of the MOST mission ( 1/101.4 minutes ). The other three fre- 
quencies correspond to the values of for (n = 2, m = 3), 
(n - 1, m = 2) and (n = 1, m - 3) respectively, taking 165 jjHz 
as A in Equation |4j 

The same analysis was performed for the photometric time se- 
ries. The result, shown in the right-hand panel of Figure [6l 
shares the dominant pattern with the 2D collapsed echelle of 
the window function, with ridges at the satellite orbital fre- 
quency and for the frequency spacing defined by Equation |H 
Figure [7] shows enlarged views of the panels in Figure [6l high- 



lighting a small frequency range around 55 jjHz, which is very 
close to 1/3 of the orbital frequency of MOST. A simple echelle 
diagram or comb analysis may yield some signal at this fre- 
quency spacing which is not due to asymptotic p modes in 
Procyon but rather the intrinsic modulation of stray light back- 
ground in the photometry. 

4. Conclusion 

From the two approaches presented here, one can draw two 
conclusions. First, the Procyon MOST light curve obtained in 
2004 shows some changes in the spectral power content after 
about seven days of observation (which were the first seven 
days of observation of the MOST mission) that covers the fre- 
quency range relevant to the p-mode search and thus requires 
special care in this search. A counterpart to this change in the 
CCD temperature telemetry indicates that it is not intrinsic to 
Procyon. 

Second, the echelle-diagramme technique must be applied to 
Procyon MOST data with care. Despite an excellent observa- 
tion window, the orbital artifacts in the data introduce appar- 
ent spacings close to the p-mode spacing expected for Procyon 
even though the gaps are very short and their signature in the 
window spectrum is very weak. However, they are unexpect- 
edly not negligible when searching for regularly-spaced peaks 
in the Fourier spectrum, reinforcing the need for a careful anal- 
ysis when looking for regular spacings in the signal spectrum 
even in the case of a very good observation window. 
In conclusion, we do not find any evidence for a stellar oscilla- 
tion signal in the MOST photometry, and the issue of the nature 
of luminosity oscillations in Procyon remains open. 
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Fig. 6. Left: Collapsed echelle diagramme for the Procyon window function. It is dominated by a peak at ^ 165 //Hz corresponding 
to the orbital frequency, and other peaks following Equation^ Right: Same diagramme for Procyon data, with similar patterns 
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Fig. 7. Respective enlargements of the Figure [6] The discontinuous appearance is due do the discrete frequency sampling of the 
spectra. 
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